The hysteresis loop area of the Ising model 
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The hysteresis of the Ising model in a spatially homogeneous AC field is studied using both 
mean-field calculations and two-dimensional Monte Carlo simulations. The frequency dispersion 
• and the temperature dependence of the hysteresis loop area are studied in relation to the dynamic 

' symmetry loss. The dynamic mechanisms may be different when the hysteresis loops are symmetric 

or asymmetric, and they can lead to a double-peak frequency dispersion and qualitatively different 
temperature dependence. 
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When a cooperative many-body system, such as a magnet, is placed in an oscillating external perturbation (such 
as a magnetic field), it may show also oscillating dynamic response. This response usually lags in time, creating a 
^ ' hysteresis loop with nonzero area. This phenomenon exists widely in, e.g., magnetic systems and ferroelectric systems 
c/2 I [1] , and has been arousing great interest for its important technical application and intriguing physics [2-4] . 

The recent theoretical [5-13] and experimental [14] studies on the hysteresis (also see Ref. [4] and references therein) 
focus on two topics: the dynamic symmetry breaking, and the area of the hysteresis loop. The first phenomenon is due 
to the competing time scales in such nonequilibrium systems [4] : The hysteresis loop loses its symmetry when the time 
" ■ period of the oscillating external perturbation becomes much smaller than the typical relaxation time of the system. 
d On the other hand, the interesting variance of the hysteresis loop area with such parameters as temperature and 
^ oscillation frequency can also be attributed to the time scale competition. For example, in the frequency dispersion of 
I \ the loop area of the Ising model, the frequency ujq/2'k giving the maximal area corresponds roughly to the point where 
' O ■ a resonance occurs. When an Ising system is placed in an AC field, the dynamics may consist of domain nucleation 
^ \ and/or domain growth [3]. The nucleation rate of new domains can be predicted by a characteristic time r„, and the 
domain growth rate is also linked with a characteristic time Tg. The resonance occurs when the time period of the 
external perturbation is comparable to either one of these time scales or a combination of them. As is shown below, 
the details of this dynamic time scale competition necessarily rely on the dynamic phase of the system. 
' ' ■ In the present work, we hope to help clarify the relationship of the two mentioned-above topics in the framework 
^ ' of the Ising model, with mean-field (MF) calculations and two-dimensional Monte-Carlo (MC) simulations. (The 
frequency range that receives the most attention here is within the discussion of the previous works using the same 
\ methods [4].) When the loops are symmetric, the system dynamics is controlled by a domain nucleation-and-growth 
OO • mechanism. When the loops are asymmetric (especially when the magnetization is well above or below zero) , through- 
' out the system evolution we can observe most spins being in the same direction. The dynamics of the remaining spins 
\l in the opposite direction may be described mainly by the domain nucleation mechanism. In the following, we shall see 
that the variance of the loop area with frequency, field amplitude and temperature strongly depends on the dynamic 
mechanism, which is determined by the loop symmetry. 

The model. Before the results are presented, we first describe the model. (1) The evolution of the magnetization 
M in the mean-field Ising model is determined by the following equation [4], 

S ■ dM ^ fM + H(t)\ 

o 

^ (2) In the MC simulation, the Hamiltonian of the two-dimensional Ising model in a spatially homogeneous field H (t) 
can be written as 

Hi{a^}■,t)^-JJ2<^^<^J-Hit)J2<^^■ (2) 

The magnetization is obtained as M = (ai). The MC simulation goes as the following: The coupling constant J 
and the Boltzmann constant kg are both taken as 1. On a two-dimensional iV x (in the present work, 100 x 100) 
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lattice with periodic boundary condition, at each time step a spin aj is randomly chosen and the probability that it 
is flipped is [15] 



W [cTi ^ ai) = ^ exp 



-—n{{ai,ajjii};t) 



(3) 



where ai = ±ai and Q is the normalization factor. Then in the next time step the field is updated and another spin 
is picked at random. A MC step consists of N x N such unit procedures. The unit time is chosen to be one MC step, 
with the time resolution being l/N"^. The initial state is always 80% randomly chosen spins up (MC simulation) or 
M = 1 (MF calculations). In both studies, the system evolves into equilibrium after a zero-field relaxation. Then, 
an AC field, H (t) = Hq sin cot, is applied. The measurement of the hysteresis loop always begins after a number of 

introductory cycles, and the symmetry of the loop can be characterized by the order parameter Q = Jq'^^'^ Mdt. 

The phase diagram of the Ising model in an oscillating field has been extensively studied [4,12,16,17]. With regards 
to the area scaling of the hysteresis loops of the Ising model, there are already MC and MF results, and a detailed 
review can be foimd in Rcf. [4]. Generally speaking, the area can be written as ^4 = Asta + Akim where Agta is 
the possibly nonzero static contribution existing even in the quasistatic limit, and Aj^in is the kinetic contribution 
assuming the form of 



HST-0g[u{w,Ho,T)]. (4) 



Here g (uj) has been believed to be a single-peak function, and u> (w, Hq,T) has been supposed to have the form of 
uj/ {H^I^T^^} There have been extensive theoretical efforts to determine the exact form of giuj), and the values and 
the physical meaning of the exponents. As we shall see below, the dynamic phase transition may lead to a double-peak 
frequency dispersion and a piecewise analytic function of temperature dependence. 

The frequency dispersion of the hysteresis loop area. Fig. 1 (a) shows a typical result of the two-dimensional MC 
simulations, at temperature T = 1.0 < Tc- (When T <Tc, a, dynamic symmetry loss can be observed as Hq decreases 
from 4 to 0. [4]) As is clearly shown in the curve with Hq <2, two distinct peaks can be observed. We name the left 
one as peak I and the right one as peak II. It is found that peak I occiirs in the range where the hysteresis loops are 
symmetric and peak II is in the range of asymmetric loops. As the frequency is increased from peak I to peak II, a 
dynamic symmetry loss occurs. 

Fig. 1 (b) shows a typical result of the MF calculations, at temperature T = 0.5 < Tc. Similarly, a double-peak 
function can be observed. However, there is a very important difference. Given relatively small values of Hq, it is 
possible that there is only peak II in the MF calculation. This is because, the equilibrium magnetization can be 
obtained by solving the following equation, 

-M + tanhf^^^^) =0. (5) 



With H small enough and T < Tc = 1, there can be two stable solutions to Eq. (5), corresponding to two values of 
stable equilibrium magnetization (a positive one and a negative one). This means that the hysteresis loops can be 
asymmetric even in the quasistatic limit. By contrast, according to Ref. [4], in MC simulation the hysteresis loop is 
always symmetric in the quasistatic limit, as a result of fluctuation. 

In the above discussions we provide evidence of two peaks I and II, which correspond to the resonance of symmetric 
and asymmetric hysteresis loops respectively. In the following we give a general explanation of the physics origin of 
this observation. The existence of two peaks clearly indicate two time scales, corresponding to two different dynamic 
mechanisms. When the loops are symmetric, both the initial domain nucleation and the late stage domain growth are 
at work. As suggested by Liu et a,l. [8,9], the observation of peak I means that a third time scale ti can be defined as a 
combination of Tn and Tg, and the resonance occurs as 27r/a;i ~ tj. Since for the asymmetric loops the magnetization 
can be always well above or below zero, at any time during the system evolution we can observe most spins having 
the same direction. The late stage domain growth is relatively inhibited, and thus, the time scale th corresponding 
to peak II shall be mainly determined by t„. 

Compared with peak II, the resonance at peak I has been relatively well studied in the previous works. Some 
illustrations of the resonance at peak I can be found in, e.g., Ref. [11,9]. In the following we focus on peak II, as 



^An exception is that, in the MF Ising model, with temperature T above Tc = 1, and Ho very small, it is analytically obtained 
A ~ HlT-'-uj/ (g^ + a;2) , where e = (T - 1) /T. [4] 
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illustrated in Fig. 2 (a) and (b). Actually, the existence of peak II can be easily understood in the MF Ising model. 
As is mentioned above, when H is small enough, there can be two stable solutions of Eq. (5) . Thus, with oj ^ 0, 
the hysteresis loop reduces to a curve (the stable solution of Eq. (5), as shown in Fig. 2 (b). In the other limit, with 
w — »■ 00, the system cannot respond to the external field, and the hysteresis loop becomes a horizontal line. Thus, 
it is straightforward to predict the existence of a peak in the intermediate region, where the phase lag between the 
external field and the system response creates a hysteresis loop with nonzero area. Note that the observation of peak 
II requires temperatures lower than the static critical point and small enough values of Hq. When Hq <^T ^ 1, we 
can give a quantitative description of the system behavior by solving the MF equation (1) analytically. We suppose 
1 — M — > 0, and obtain 



In the following we report some important differences and similarities of peak I and II, as summarized from Eq. 
(6) and the numerical results in Fig. 1. Differences: (1) The temperature dependence of the height of peak II in the 
MF Ising model is obtained in Eq. (6), and is different from the previous T~^/^ prediction of peak I [17]. (2) In the 
MF Ising model, the maximal area ^njax '^.t peak II grows with Hq as Hq, while the maximal area A]^^^^ at peak I 
has been predicted to grow with Hq linearly in the previous studies [17]. This difference is explained by observing the 
variance of the loop shape with Hq, as illustrated in Fig. 3(a) and (b): At peak I, only the width of the loop increases 
with Ho (that is why Aj^g^^ grows linearly with Hq), while at peak II, the loop is expanding in two directions with 
^max growing as H^. 

Similarities: (1) In both MF calculations and MC simulations, it is found that, as w — > oo, the area decays as lo~^. 
(With respect to peak I, this is in accordance with the previous MF result and the work of Rao et al. [18,19] on the 
(<I>^)^ and ($^)^ model, but not the previous MC result of the Ising model [4,17], which indicates an exponentially 
decaying function of g [u (w. Ho, T)] in Eq. (4).) (2) Independent of Hq and T, peak II is always observed to be at 
(or very close to) cuu = 1. With regards to peak I, in both MC simulations and MF calculations we find that as , 
also approaches 1, which is different from the previous predictions. According to Ref. [4,17], the uji will also tend to 
infinity as Hq — > oo. But it is not what we observe in Fig. 1, which shows that, as the field already far exceeds the 
spin- spin interaction, the time scale of the system is no longer sensitive to the value of Ho- 

In the following we turn to study the temperature dependence of the loop area with fixed field amplitude and 
frequency [17,13]. Here our motivation is quite similar to that of the above discussions of the frequency dispersion. 
When Ho < 4 (MC) or Hq < 1 (MF), a dynamic symmetry loss can be observed as T decreases (as can be predicted 
from the phase diagram previously obtained [4,12,16,17]). Thus, a simple scaling function is not likely to exist for the 
loop area, since there are different dynamic mechanisms of the symmetric and asymmetric loops, and different time 
scale competition. This is supported by the MF and MC results. 

A typical MF result is shown in Fig. 4 (a). The order parameter \Q\ > for lower temperature and [Q] = 
for higher temperature, and at the dynamic critical point, a dynamic symmetry loss occurs. In Ref. [13], using the 
same methods as the present work, Acharyya found that the area becomes maximum above the dynamic transition 
point. Here, it is clear that the temperature dependence of the loop area assumes different functions for the [Ql > 
states and [Ql = states. These different functions are separated by the dynamic critical point and the first order 
derivative, dA/dT, is not continuous at the dynamic critical point. Thus, the temperature dependence is a piecewise 
analytic function. A typical MC result is shown in Fig. 4 (b), and it is roughly similar to the MF result. Although 
we do not observe a notable discontinuity of the first order derivative, dA/dT, it is very obvious that the second 
order derivative, d^A/dT"^, changes sign at the dynamic critical point. When iJo > 4 (MC) or TJq > 1 (MF), the field 
amplitude exceeds the spin-spin interaction, and the hysteresis loops are always symmetric [4,12,16,17]. In this case, 
we observe that the loop area decreases monotonically as the temperature grows. 

To summarize, the hysteresis of the Ising model in an AC field, H (t) = Hq sin [ujt) is studied using both mean-field 
calculations and two-dimensional Monte Carlo simulations. The frequency dispersion and the temperature dependence 
of the loop area are studied in relation to the dynamic symmetry loss. The dynamic mechanisms are different when the 
hysteresis loops are symmetric or asymmetric. For symmetric loops, the dynamics is a combined domain nucleation- 
and-growth process. By contrast, for asymmetric loops well above or below the M = line, the dynamics may 




This equation can be exactly solved, and the loop area is obtained as 



A=- 




(6) 
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be mainly domain nuclcation. This framework is part of basic current knowledge of hysteresis phenomena, and the 
observed frequency and temperature dependence of the loop area is consistent with it. Double peaks can be observed 
in the frequency dispersion, and the temperature dependence is possibly a piecewise analytic function. Interestingly, 
the shift of the dynamic mechanism with the symmetry loss is also found in the mean-field calculation, and some 
striking similarities are observed (for example, the same position of peak II). Although the present work deals with a 
model spin system, the topics discussed have general meaning. Surely some quantitative details like the position of 
the peaks rely on the model setting, but we believe that the physics of the conclusions is not limited to the specific 
system studied here, and can be predicted for more general systems. 

We thank Hao Yu for helpful discussions. This work is supported by the Natural Science Foundation of China 
(50332 020, 10021001) and National Key Projects for Basic Research of China (2002CB61 3303). 
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Fig.l The frequency dispersion of the loop area A [lo). (a) The MC results with AC field at temperature T = 1.0; 
(b) The MF results with AC field at T = 0.5. The dashed line, corresponding to io~^ , serves as a guide to eye. 

Fig. 2 The illustrations of peak II with (a) MC simulations at T = 1.0 and (b) MF calculations at T = 0.5. 

Fig. 3 (a) Some typical hysteresis loops at peak I, with i?o = 3,6, ...,30 from the innermost to the outermost loop, 
(b) Some typical hysteresis loops at peak II, with Hq — 0.05, 0.15, 0.25. All the results are obtained from MF 
calculations at T = 0.5. 

Fig. 4 The temperature dependence of the loop area A and the order parameter Q obtained from (a) MF calculations 
and (b) MC calculations, with the AC field. 
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